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Abstract 

In the present paper, the stability of Coalescence Hidden variable Fractal Interpo- 
lation Surfaces(CHFIS) is established. The estimates on error in approximation of the 
data generating function by CHFIS are found when there is a perturbation in inde- 
pendent, dependent and hidden variables. It is proved that any small perturbation in 
any of the variables of generalized interpolation data results in only small perturbation 
of CHFIS. Our results are likely to be useful in investigations of texture of surfaces 
arising from the simulation of surfaces of rocks, sea surfaces, clouds and similar natural 
objects wherein the generating function depends on more than one variable. 
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1 Introduction 



The theory of fractal interpolation has become a powerful tool in applied science and engi- 
neering since Barnsley pQ introduced Fractal Interpolation Function (FIF) using the theory 
of Iterated Function System (IFS). Massopust [8] extended this concept to Fractal Interpo- 
lation Surface (FIS) using IFS wherein he assumed the surface as triangular simplex and 
interpolation points on the boundary to be co-planar. In view of lack of flexibility in his 
construction, Geronimo and Hardin [6] generalized the construction of FIS by allowing more 
general boundary data. Subsequently, Xie and Sun [10] used bivariate functions on rectan- 
gular grids with arbitrary contraction factors and without any condition on boundary points 
to construct Bivariate FIS. Dalla [5] improvised this construction by using collinear bound- 
ary points and proved that the attractor is continuous FIS. However, all the constructions 
mentioned above lead to self-similar attractors. 

A non-diagonal IFS that generates both self-affme and non-self-affme FIS simultaneously 
depending on the free variables and constrained variables on a general set of interpolation 
data is constructed in [2]. The attractor of such an IFS is called Coalescence Hidden- 
variable Fractal Interpolation Surface (CHFIS). Since the CHFIS passes through the given 
data points, any small perturbation in the data points results in the perturbation of the 
corresponding CHFIS. 

The construction of a Coalescence Hidden-variable Fractal Interpolation Function (CHFIF) 
of one variable and investigation of its stability is studied in 0, H] . A CHFIF is an important 
tool in the study of highly uneven curves like fractures in rocks, seismic fracture, lightening, 
ECG, etc. However, it can not be applied for the study of highly uneven surfaces such 
as surfaces of rocks [10], sea surfaces [9], clouds [11] and many other naturally occuring 
objects for which the generating function depends on more than one variable. A CHFIS is 
a preferred choice for the study of these naturally occurring objects. The quantification of 
smoothness of such surfaces in terms of Lipschitz exponent of its corresponding CHFIS is 
investigated recently in [7]. The purpose of the present paper is to investigate the stability 
of such CHFIS. The estimates on error in approximation of the data generating function by 
CHFIS are found individually when there is a perturbation in independent, dependent or 
hidden variable. These estimates together give the total error estimate on CHFIS when there 
is perturbation in all these variables simultaneously. It is proved that any small perturbation 
in any of the variables of generalized interpolation data results in only small perturbation 
of CHFIS. Unlike the case of CHFIF, the stability of CHFIS is studied with respect to 
Manhattan metric and requires the perturbations in generalized interpolation data to be 
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governed by an invariance of ratio condition. Our results are likely to find applications in 
texture of surfaces of naturally occurring objects like surfaces of rocks, sea surfaces, clouds, 
etc.. 

A brief introduction of CHFIS is given in Section [2j In Section [3j some auxiliary results 
that are needed to establish the stability of CHFIS are derived. Our main stability result 
is established in Section H] via stability results found individually for perturbation in inde- 
pendent variables, the dependent variable and the hidden variable. Finally, these results are 
illustrated in Section [5] through simulation of a sample surface with given data as well as 
with perturbed data obtained by small variations of concerned variables in the given data. 



2 Preliminaries 

For the given interpolation data {(xj, yj, Zij) : % — 0, 1, . . . , N and j — 0,1, ... , M}, where 
Zij G (a, b) and — oo < a < b < oo, consider the generalized interpolation data A = 
{(xj, Dj, Zij, tij) : i = 0,1, ... ,N and j = 0,1, ... , M}, where tij G (c, d) and — oo < c < 
d < oo. Set I = [x , x N ], J = [y , y M ], S = I x J, D = (a, b) x (c, d), I n = [sc n _i, x n ], J m = 
lUm-i, Um] and S ntTn = I n x J m for n = 1, . . . , N and m = 1, . . . , M. Let, the mappings 
4> n : I — > I n , ip m : J — > J m and F n ^ m : S x D — > D for n = 1, . . . , N, and m — 1, . . . , M be 
defined as follows: 

4> n {x) = X n _i + — -2=i (x - x ), 

i I \ , Dm Dm— I , \ 

^m{y) = y m -i + [y - yo), 

Vm - yo 

-^n,m(x, y, Z, t) ifin.m % fn,m V ^n.m z Pn,m t ~\~ Qn,m ^n,mi 

Cn,m X -|- fn,m y ^n^m t ~\~ Qn,m %V ^n,m)- 

Here, a n , jm and ^ n ,m are free variables chosen such that \a n , m \ < 1 and |7 n ,m| < 1- \P n ,m\ is 
a constrained variable chosen such that \(3 n ,m\ + \ln,m\ < 1- Let the function F n ^ m satisfy the 
following join-up condition: 



— l,m— 1 ) 
,m— 1) tn,m— 1 ) 

F n ,m\ x 0i VM, Z0,Mj to,M) — (Vl,m)in-l,m) 
Fn,m( x N-,UMi z N,Mi^N,m) = ( z n,m,tn,m] 



(2.1; 
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Using the condition (J2HD, the values of e n>m , f n>m , g n>m , k n>m , e n>m , f n>m , g n>m and k n , m are 
determined as follows: 



9n,m 
fn,m 




— l,?n- 


1 &n — 1,171 2?i,m — 1 -\~Zn,m Otn,mZeva Pn 7 mteva 






— l,m- 


(x -x N )(y -y M ) 

1~ ^n.m-l— an,m(zo,0 — z N,o)~ 0n,m(to,Q— *JV,o) — Sn.m^O - 


-xn)vo 


Zn 


— l.m- 


(X0-Xiv) 

1 — z n-l,m— a n , m (2o,u— 2q,m)— /3n,??i(io,0— to,Af) — 3?i,m(y0 


-Vm)xo 






(vo-Vm) 




h 

">n,m 




m 


£n,m%N ~ fn,mVM ~ (%n,m z N,M ~ Pn,mtN,M 


- g n ,mXNyM 


9n,m 
&n,m 
fn,m 


tn- 


-l,m— 


1 Cn — 1 ,m tn,m — 1 ~\~in } m Tn,m ^eva 




tn- 


- l.m — 


(x -x N )(y -yM) 
1 — *n,m-l — 7n,m(*0,0— *iV,o) — Sn.mC^O— Xff)yo 




tn- 


- l.m — 


(x -x N ) 

1 — tn-l,m— 7»,m(to,0— *0,Af) — §n,m(j/0— Vm)^>0 








(vo-Vm) 




k 


tn 


m 


&n,mXN — fn,mUM ~ ln,m~tN,M ~ Qn^NVM 





} (2-2) 



where, z eva = zn,m — Zjv,o — ^o,a/ + ^0,0 and t eva = tjf,M ~ *at,o — ^o,M + £0,0 • 
Now define the functions G n ^ m {x, y, z, t) and u nim (x, y, z, t) as 

{F n+1>m (xo,y,z,t), x = xn, n = 1, . . . , N - 1, m = l,...,M 
^k,m+i(a:,2/o,2,*)> V = VN, n = l,...,N m = l,...,M -1 
F n>m (x,y,z,t), otherwise. 

and u n>m (x,y,z,t) = (cf) n (x),^ m (y),G n>m (x,y, z,t)). Then, {M 4 ,u; nim : n = l,...iV; m = 
1, . . . M} constitutes an IFS for the generalized interpolation data A. It is known [2] that 
there exists a metric r on M 4 , equivalent to the Euclidean metric, such that the IFS is 
hyperbolic with respect to r and there exists a unique non-empty compact set G C M 4 

N M 

such that G = (J (J cj n)m (G). The set G is called the attractor of the IFS for the given 

n=l m=l 

interpolation data. Further, G is the graph of a continuous function F : S — > M 2 such 
that = (zij, tij) for i = 0, 1, . . . , N and j = 0,1, ... ,M i.e. G = {(x, y, F(x, y)) : 

(x,y) G S and F(x,y) = (z(x,y),t(x,y))}. 

Definition 2.1 Let F{x,y) be written component-wise as F(x,y) = (Fi(x,y), F 2 (x,y)). The 
Coalescence Hidden-variable Fractal Interpolation Surface ( CHFIS) for the 

given interpolation data {(xi, yj, Zij) : i = 0, 1, . . . , N and j = 0, 1, . . . , M} is defined as the 
surface z = F\(x, y) in R 3 . 

It is easily seen that the function F(x, y) described above satisfies 

F{x,y) = G rhrn { ( l ) ; i \x),^\y),F^-\x),^\y))) 

for all (x,y) G S n>m , n = 1,2, . . . , N and m = 1,2, . . . , M. Consequently, if the func- 
tion G n , m (x,y,z,t) is written component-wise as G n , m (x,y, z,t) = (G} im (x,y),G 2 nm (x,y)), 
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then 



Fx((f) n (x),ip m (y)) = G 1 ntm (x 1 y 1 F 1 (x,y),F 2 (x,y)) = a n>m Fi(x,y) + n ,m F 2 (x,y) +p n ,m(x,y) 
F2(<Pn(x),ip rn (y)) =G% m (x,y,F 1 (x,y),F 2 (x,y)) = 7„, m F 2 (x, y) + q n , m (x, y) 



where, 



Pn,m(.X, y) — 6n,m 2< ~i~ fn,m 2/ ~i~ 9n,m Xy -\- k n m 
Qn,m(.X, y) — Cn,m X + fn,m y ~\~ 9n,m Xy + k njm . 



(2.3) 



3 Some Auxiliary Results 

In this section, we develop some results that are needed in the sequel for investigating the 
stability of CHFIS in Section H 

Let — oo < Xq < x\ < . . . < x* N and — oo < y$ < y\ < . . . < y* N be a bounded set of 
real numbers in x and y axis. Denote I* = [xq,x* n ], J* = [yl,y* M ], S* = I* x J* I* = 
[<_!, <], J* m = [y* m ^ y* m \ and S* >m = I* x J* m for n = 1, . . . , N and m = 1, . . . , M.. Similar 
to 4> n and ip m in section [2J we construct homeomorphisms 0* and ^)* m where 0* : I* — »■ I* 
and : J* — > J^. Now, define the map R n .m '■ S n ,m — > S* m by 



t-> I \ I * i n ~n— 1 / \ * . am om—L / \ 

Rn,m[X, y) = I x n _i H (x - X„-i), H (y - J/ m _i) 

\ *^n— 1 l/m Vm—X 

and the linear pi6CGwis6 nicip 7? ; — *S* £is 

R(x,y) = R n , m (x,y) for all x E I n and y e J m (3.1) 

Similarly, the maps K n ^ m : S 1 * m — ► S n ,m an d K : S* ^ S are defined as 

#n, m (x , y ) = I Xn-l + „ - Z„-i) , Vm-1 + \V ~ Vm-l) 

and 

K{x*,y*)=K n , m (x*,y*) for all x* G I* and G J,* n . (3.2) 

It is easily seen that K = R~ l . Set £ n , m (x, y) = {(j) n {x), if> m (y)) and £* im (z, y) = (^*(z), V m {y))- 
We assume the following invariance of ratio condition for any two sets A = {(cCj, j/j, Zjj, t^j) : 
i = 0, 1, . . . , N and j = 0, 1, ... , M} and A* = {(x*,y*, z itj , t it j) : i = 0,1,..., N and j = 
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0, 1, . . . , M} of the generalized interpolation data points 



(x - X N ) _ - X n ) 



and 



(yo ~ Vm) _ (y, 




(3.3) 



x X JV/ l x n-l x n) 



(vo - v* M ) (y*„ 



By ( 13. 3p . we observe that for n = 1, 2, . . . , N and m = 1, 2, . . . , M, 



(Ro^ m oK) {x*,y*) 
F„, m (K(x*,y*),z,t) 



(3.4) 



Thus, the dynamical systems {S; £ n ,m} and {S*; £* m } are equivalent. Using this equivalence 
of dynamical systems, we first prove the following proposition needed for establishing the 
smoothness of CHFIS in Proposition 13.21 

Proposition 3.1 Let A = {(xi,yj, Zij,tij) : i = 0,1,..., N and j = 0,1,..., M} and 



alized interpolation data for the IFS {3ft 4 ; w n>m (x, y,z,t),n = 1, . . . , N; m = 1, . . . , M} and 



free variables and constrained variable. Let the points of A and A* satisfy the invariance 
of ratio condition (13.31) . Then, F, as defined in Section [U is the CHFIS associated with 
{3ft 4 ; u) n>m ,(x, y,z,t), n — 1, . . . , N; m — 1, . . . , M} iff F o R' 1 is the CHFIS associated with 
{$¥;LO* m (x*,y*,z,t), n = l,...,N;m = 1,...,M}, where R is defined by (EOT) . 

Proof Let F(x, y) be the CHFIS associated with the IFS {9ft 4 ; LO n>m (x, y, z, t), n — 1, . . . , N, 
m = 1, . . . , M}. It follows from (El) that, 



K >m (C;^*,y*), FoR-\C n ^(x\yl))=F nim (KoC n ^(x*, y *), FoR~\C n ^(x\yl)) 



Thus, FoR- 1 is the CHFIS associated with {3ft 4 ; uj* n m {x* ,y* , z, t), n = 1, . . . , N, m = 1, . . . , M}. 
Conversely, assume that the above identity holds for CHFIS F o R~ l . 



A* = {(xly^Zij^ij) : i = 0,1,... 



N and j = 0,1,..., M} be any two sets of gener- 



{$l' l ;ulm(x*,y*,z,t),n = l,...,N;m 



M} respectively, with the same choice of 



F n , m (e~i o K(x*,y*),F{d o K(x*,y*))) 
F^^ix), ^-\y), F{4>-\x), ^-\y))) 
F(x,y) 

FoR-\ x *,y*) 
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Then, for x G I n and y G J m , 

F n,m (CK^. V), F (£nJm( X , V))) = F n,m ° C^O, V)i ^(CmO, 2/))) 

= (C- ° R{*> v)' Fo R~\C n ^ O R(X, y ))) 

= FoR~\ x \y*) 
= F(x,y) 

Hence, F(x, y) is the CHFIS associated with the IFS {!R 4 ; uj n>m {x, y,z,t), n = 1, . . . , N, and 
m = l,...,M}. ■ 

2 

Let X = (xi,X2) , Y = (2/1,2/2) ^ I^ 2 and cIm{X, Y) = E l x « — Vi\ be the Manhattan metric 

i=l 

on I 2 . A function F : M 2 — > M is said to be a Lipschitz function of order 5 (written as 
F G Lip 5) if it satisfies the condition \F(X) - F(Y)\ < c [d M (X, Y)] 5 , 5 G (0, 1] and c is a 
real number. The real number 5 is called the Lipschitz exponent. 

Define 

Qn(F 1 (X)) := E E 



where, Yq (X) 

5 /YOri,...,r n ,si,...,s n \ J 



Yj E ( XSr 1 ,...,r n ,s 1 ,..., Sn (X) \Sr^,.'.,rn,n,-';Z\ 

i,...,r»=l si,...,s re =l \ 



1, X G 5Vi,...,r n ,si,...,s n 

0, otherwise. 



and 6 ri ,..., r „, fll ,..., Sn = / / Fi(x,y) dy dx. 



Iri r„ ^sj,...,sn 



It is known [7] that Qn(-FipT)) converges to Fi(X) uniformly with respect to Manhattan 
metric. Using this fact and finding a bound on maximum distance between Q n (Fi(X)) and 
Q n (Fi(X)), the Lipschitz exponent of CHFIS F x is found in the following proposition when 
p n:m and q n>m , given by ([23D, belong to Lip 1. 

Proposition 3.2 Lei Fi(x,y), < 2,2/ < |, 6e £/ie CHFIS for the interpolation data A = 
{(xi,2/j,^j,tij) : i,j = 0,1,..., iV} w/iere, x = 2/0 = and x N = y N = 1/2. T/ien, 
Fi G Lip 5 for some 5 G (0, 1]. 

Proof Using the values of e n , m , / n , m , c7n, m , fc n>m , e n , m , / n , m , g n , m and fc n , m as defined 
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in (12.21) . we find that, for < x < x < \ and < y < y < \ , 

\p n ,m(x,y) -p n ,m(x,y)\ < 8 • [(1 + a) Z max + (3 T max ] ■ [\x - x\ + \y - y\] 
\q n , m (x,y) - q nim (x,y)\ < 8 • [(1 + 7) T max ] ■ [\x -x\ + \y- y\] 

where, Z max = max{|^j - z k ,i\ } and T max = max{|t iji - t k>l \ } for i,j,k,l = 0,...,N. 
Therefore, the bound on \Q n (Fi(X)) — Q n (Fi(X))\ is obtained as 

|Q n (F 1 (X))-g„(F 1 (X))| 

< M[d M (x,x)]{[i + (&) + ... + ( n)( m - 2 )] 
+ \p\ ■ [( r)[i + . . . + ( n)( m - 4 >] + ( f ) 2 [i + ... + ( nY m ~V] + ... + ( f r~ 3 ] j (3.5) 

The above inequality (13. 5p is similar to the Inequality 4.1 of [7] except that, in this case 
we have, Q = Q = G and T = T. By employing the same technique of proof here and 
using the inequalities < — (|a;| + \y\) T log(|ar| + \y\) < ^ and < (|x| + |y|) T (log(|x| + 
\y\)) 2 < 7^2? we find that, for different magnitudes of and f (0 > 1, = 1 or < 
1; f > l,f = 1 or f < 1), the bound on IQ^F^X)) - Q^F^X))] is obtained as 
\Q n (F 1 (X) - Q n {F 1 {X))\< M[d M (X,X)] 5 for some 5 E (0,1]. Using the fact Q n (Fx (X)) 
converges to Fi(X) uniformly with respect to Manhattan metric, we get Fi G Lip 5. ■ 

Corollary 3.1 Let F\ be the CHFIS for the generalized interpolation data A = {(xi, yj, Zij, tij) : 
i, j = 0, 1, . . . , N} . Then, there exist constants K and 5, independent of z n , m and t n>m , such 
that 

\Q n (F 1 (X)) - Q n (Fi(X))\ < MK [d M (X,X)] S (3.6) 
where, X = (x, y) and X = (x, y). 

Proof Define a linear homeomorphism R : S — ► [0, ~] x [0, |] that transforms the given 
interpolation data A to the data Aq = {(0, 0, z ,0j t ), (a;*, 0, Zifl, £i,o), • • • > (\, 0, zn,o, ^,0); 
(0,y*,zo tl ,t 0>1 ),...,(0, l,zo tN ,t 0)N ), ...,(§, l,z NiN ,t N)N )}. Proposition (jH^J) applied on the 
data A* gives \Q n {Fi o R~\X*)) - Q n (iq o R^(X*))\ < M[d M {X*,X*)] S . Now, using 
Proposition ( 13. ip with A and Aq, the corollary follows. ■ 
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4 Stability of CHFIS 



To prove the main stability result of CHFIS, we need to investigate its stability with respect 
to perturbations in independent variables, the dependent variable and the hidden variable 
of the generalized interpolation data. We develop these results first. 

The following theorem gives the effect on stability of CHFIS due to perturbation in inde- 
pendent variables: 

Theorem 4.1 Let F\ and G\, with the same choice of free variables and constrained variable, 
be CHFIS respectively for the generalized interpolation data A = {{xuVj, Zi j,ti j) : i,j = 
0, 1, . . . , N} and A* = {(x*, y*, Zij, Uj) : i, j = 0, 1, . . . , N} that satisfy the the invariance of 
ratio condition (13.31) and S* C S, where S* = [xq,x* n ] x [y$, y* M ] and S = [xq,Xn] x [yo,yM]- 
Then, 

ll^i-GiiL 



< M 



Wl + (1 + ot) 



(l-a)(l- 7 ) {I -a) 



max \ (\x n - x* n \ + \y m - y* m \) S : n, m = 0, . . . , N 



(4.1; 



Proof Let the function G(x, y) corresponding to the generalized interpolation points A* 
be defined as G(x,y) = G* ffn (^(x), ^(y), ^{y))) for x E [<_!,<] and 

V e [y* m -i,V*J- By (O and (JH, we observe, 

d (R(Mx), Mv))) = Km (x, y, g ° r(x, y )) 

= a n ,m (Gi o R)(x, y) + f3 n , m (G 2 o R)(x, y) 
+ e njm x + fnjn y + g n , m xy + k 

n,m 

G 2 (R((f)n{x),^ m (y)) = Fn.rn (», D, G O R(x, y)) 

(G 2 oR)(x,y) + x + f n ,m y + g n,m Xy + kn^m 

Thus, we have, 

|Gi(W„(x),^ m (y))) - F 1 ( ( f )n (x)^ m (y))\ 

< \a n ,m\ \G\ o R(x,y) - F x {x,y)\ + |/3„, m ||G 2 o - F 2 (x, y)| 

< |a| |Gioi2(x,y)-Fioi2(x,y)| + |a| |F X o R(x, y) - F t (x, y)\ 
+ \/3\ \G 2 oR(x,y)-F 2 (x,y)\ 

Using Proposition 13.21 and Corollary 13.11 it follows that there exist constants M and 5 
which are independent of z n)m and t n>m such that \Fi(X) — F\(X)\ < M dM(X,X) 5 where, 
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X = (x, y) and X = (x, y). Since (x, y) G S implies (x, y) G S n>m for some n,m — 1, . . . , N, 
it is easily seen that R(x,y) G S* m which in turn gives cIm(X, R(X)) < max{(|a; n — x* n \ + 
\Vm — Vm\) '■ n,m = 0,1, ... , N} . Thus, the last inequality reduces to 
\\G l oR-F 1 \\ OD < \a\ \\G 1 -F 1 \\ O0 + \/3\ \\G 2 oR-F 2 \\ 00 

+ \a\ Mmax.U\x n - x* n \ + \y m - y* m \) 5 : n, m = 0, 1, . . . , n\ 



Now, ||Gi - FiHoo < - F 1 o iT 1 ^ + \\Fx o R~ l - FxW^. Therefore, 

— a) 
[1 + a) 



\Gx — -Fllloo < 1~J~ \ 11^2 ° R — -^2 1 loo 

(1 - a) 



f ,l-a) 



Mmax {(\x n -x* n \ + \y m - y*J) 5 : n, m = 0, 1, . . . , iv} (4.2) 



In order to find an upper bound of \\G 2 ° R — F 2 \\oo in (14. 2p , we observe that, 
\G 2 o R((f> n (x),il) m (y)) - F 2 {(j)n{x),'4>m{y))\ 

< | 7 | I \G 2 o R(x, y)-F 2 o R(x, y) \ + \F 2 o R(x, y) - F 2 (x, y) \ 
The above inequality holds for all n, m — 1, 2 . . . N. So, 

||G 2 o R- Falloo < |7| \\G 2 - Falloo + (7! Mmax - x* n \ + \y m - y* m \) & : n, m = 0, . . . , n\ 
Hence, 

||G 2 - FalU < ^j^T^} max - <| + \y m - y*J) 5 : n, m = 0, . . . , iv} 
Therefore, 

HG^ofl-^Hoo <M-^-j max{(|x n -x;| + |y m -^|) <5 :n,m = 0,l,...,iv} (4.3) 

Substituting inequality (14. 3 p in inequality (14.21) . we get the required bounds. ■ 

The following theorem gives stability of CHFIS when there is a perturbation in the dependent 
variable. 

Theorem 4.2 Let Fx and Gx, with the same choice of free variables and constrained variable, 
be CHFIS respectively for the generalized interpolation data A = {(xi,yj, Zij,tij) : i,j = 
0,1,..., N} and A* = {( Xi , Vj , z^, t id ) : i,j — 0,1, ... , N}. Then, 

ll^i-GilL < 4 / n 1 + Q ? maxl\zn,m-z* nm \:n,m = 0,l,...,N\ (4.4) 



[1-a) 

Proof In view of Proposition 13.11 we may assume [xq,xat] x [j/o? S/jv] = [0, |] x [0, |]. Since 
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the independent variables and hidden variable are same in both the interpolation data, the 
self affine FIS are the same i.e F 2 (x, y) = G 2 (x, y)- The value of e n>m , f n ,m, 9n,m, k nym differs 
from e* m , f^ m , g* nm , Kim as the perturbation occurs in dependent variable. 

Therefore, for n,m — 1,2, ... ,N, 
\Fi(<f> n (x),ij} m (y)) - Gi((j) n (x),ip m (y))\ 
< \a\ \Ft(x,y) - G x (x,y)\ 



[l + a) 



1, , 1, 1, 

2\x--\+2\y--\+4\xy--\ + l 



max <^ \z„ m - z* m I : n, m = 0, 1, . . . , N 



Since [2|x — ~| + 2\y — ~| + 4:\xy — ~| + l] < 4, the required bounds for stability is obtained 
from the above inequality. ■ 

When the hidden variable is perturbed, both the self affine function F 2 and CHFIS F-y gets 
perturbed. Proposition 14.11 describes the stability of self affine function F 2 when the hidden 
variable is perturbed. Using this proposition, the stability of CHFIS is described in Theorem 
(14. 3 p when the hidden variable is perturbed. 

Proposition 4.1 Let F\ and G\, with the same choice of free variables and constrained 
variable, be CHFIS respectively for the generalized interpolation data A = {(xi,yj, Zij,t it j) : 
h j = 0, 1, • • N} and A* = {{x^y^z^,^) : i,j = 0,1, . . . , N}. Then, 

\\F 2 - C2IL < Tf^j max{|t„ im - t*J : n = 0, 1, . . . , N; m = 0, 1, . . . , N} (4.5) 

Proof By Proposition 13. 1[ we may assume [xo,xat] x [2/0,2/iv] = [0, |] x [0, |]. Hence, 
\F 2 ((f) n (x), ip m (y)) - G 2 (<p n (x), ip m (y))\ 

< | 7 | \F 2 {x,y)-G 2 {x,y)\ + 

x (1 +7) {max{|f n)m - t* J :n,m = 0,1,..., N}} 



1, ,1, 1, 
2|z--| + 2|y--| + 4|zy-~| + l 



The above inequality is true for all <f> n (x) and ip m {y); n,m = 1,2, ... ,N giving the required 
bounds for ||F 2 - G 2 \\ OQ . ■ 

Theorem 4.3 Let Fi and G\, with the same choice of free variables and constrained variable, 
be CHFIS respectively for the generalized interpolation data A = {(xj, yj, Zij, t^j) : i,j = 
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0,1,..., N} and A* = {{x i ,y j ,Zij,t* iJ ) : i, j = 0, 1, . . . , N}. Then, 

8 (3 f 

\\ F i ~ G i\\oo ^ 71 \7i rmax \ \t n>m - t* n m \ : n, m = 0, 1, . . . , N } (4.6) 



(l-a)(l- 7 ) 

Proof We may assume [xq,xat] x [i/o,2/iv] = [0, |] x [0, |] by Proposition 13. 1[ Thus, 



m (y)) - Gi((f) n (x), 1p m 


(y))\ 


(x,y) - G!(x,y)\ + \(: 


i\\F2{x,y)- 






2\x-\\ + 2\y-\\, 





max{|t njm - t* | : n,m = 0, 1, . . ., N} 



Using equation (l4.5l) and since [2\x — \\ + 2 1 y — \ \ + 4|xy — ^| + l] < 4, we get the desired 
bounds in (14.31) . ■ 

Theorems (I4.ip - (I4.3I) suggest the following definition of a metric on a generalized interpola- 
tion data needed to formulate our main result on stability for CHFIS. 

Definition 4.1 Let 5(A) be the set of generalized interpolation data. The metric d(A l , A 2 ) 
on the set S(A) C 3? 4 is defined as: 

d ( Al ' A2) = 71 TTi V max {lCn - £,J : n, m = 0, 1, . . . , N} 

(1 - a) (1 -7) 



4(1 + a 



-max{|4,m - 4,ml • n,m = 0,1, ... ,N} 



+ 



1 — a t 

2(3 7 (1 + a 



Mx 



(l-«)(l- 7 ) (l_ a ). 
x max{(|x^ - x 2 | + \y x m - y 2 m \) & : n, m = 0, 1, . . . , N} (4.7) 

where A m = {(x™, yj 1 , z™j, t™ ), : i,j — 0,1, ... , N} e S^A), m = 1, 2 and satisfy invariance 
of ratio condition (13.31) . 

Using Theorems (I4.ip - (I4.3I) . the main stability result for CHFIS is now obtained as follows: 

Theorem 4.4 Let Fi and G±, with the same choice of free variables and constrained variable, 
be CHFIS respectively for the generalized interpolation data A = {(xi,yj, Zij,tij) : i,j = 
0, 1, . . . , N} and A* = {(x*, y*, z*j, t*j) : i, j = 0, 1, ... , N} that satisfy the invariance of 
ratio condition (13. 3p . Then, 

HFx-GxIL^d^A*). 
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Proof Let A* = {(x u y j} z itj , t* d ) : i,j = 0,1,..., N} and A* = {(x u y j} z^, t* tj ) : i,j = 
0, 1, . . . , N} be two generalized interpolation data having CHFIS F* and F±* respectively 
with the same choice of free variables and constrained variable as for F\ and G\. By con- 
sidering the pairs (A, A*), (A^A^) and (A^A*) and applying Theorems fl4.ip - fl4.3l ) for 
these sets of data with appropriate CHFIS, it follows that, 

||*i - Gilloo < \Wi - F*|U + - F**|U + \\F? ~ GiWoo < d (A, A*) . ■ 



5 Error bounds for a Sample Surface 



Consider the sample CHFIS (c./. Fig. [I]) generated by the data in rows 1 — 3 of Table [T] with 
0.7, (3 n , m = 0.4 and 7n,m — 0.5. Perturbed values in independent variables x n and y m 
are chosen such that the invariance of ratio condition f)3.3p is satisfied. However, perturbed 
values of dependent variable z n) m and hidden variable t n{m are randomly generated. 



Figs. 2(a) and 2(b) are simulations of CHFIS generated corresponding to perturbations in 



independent variables x n and y m (c.f. Table\l\ Case (i(a)) and Case (i(b)J). The bound on 
error \\F — G||oo {c.f. Theorem \4.1\ ) for these simulated surfaces are found to be 0.0217 and 



2.1667 when max { (| 
Tabled . 



J. > , 



+ \y m ~ Vm\) '■ n ' m = 0' 1) ^ } is 0.002 and 0.2 respectively {c.f. 



Figs. 3(a) and 3(b) give simulations of CHFIS obtained by perturbing the values of dependent 
variable z n ^ m {c.f. TableUl Case (ii(a)) and Case (ii(b)n. The bound on error \\F — G\\oo 
{c.f. Theorem \4-^ for these simulations equals 0.0227 and 2.2667 when max{|z* m — z^ m \ '■ 



n, m 



0, 1,2} is 0.001 and 0.1 respectively {c.f. Table\ 



Figs. 4(a) and 4(b) demonstrate the effect of perturbations in hidden variable t n ^ m {c.f. 
Tabled Case (iii(a)) and Case (iii(b)J). The bound on error \\F — G||oo {c.f. Theorem \4-3\j 
for these simulated surfaces equals 0.0213 and 2.1333 when max{|t^ — t^ m \ : n, m — 0, 1, 2} 
is 0.001 and 0.1 respectively {c.f. Table\E>- 



Finally, Fig. 5(a) and 5(b) give the perturbed images that is simulated by simultaneously 
using perturbed independent variable, perturbed dependent variable z* nm and perturbed 
hidden variable (c.f. Table\l\ Case (i(a)) ,(ii(a))(iii(a)) and Case (i(b)),(ii(b))(iii(b)) respec- 
tively ). The computed bound on error \\F — Cr||oo {c.f Theorem \4-4\ for these perturbations 



in all the variables is found to be 0.0657 and 6.5667 when max { (| 



X v 



+ \yr, 



241) 



n,m = 0,1,2} is 0.002 and 0.2, max{|^ m — z^ m \ : n,m = 0,1,2} is 0.001 and 0.1 and 
max{|t^ m — t^ml '■ n,m = 0,1,2} is 0.001 and 0.1 respectively. 
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Table 1: Data points and their perturbations for a sample surface 



10 
H 


(x n , y m ) 


(0,0) 


(0,1) 


(0,2) 


(1,0) 


(1,1) 


(1,2) 


(2,0) 


(2,1) 


(2,2) 



0. 

« 


z n m 


0.3 


0.5 


0.6 


0.7 


0.4 


0.6 


0.8 


0.5 


0.6 


Dat 


tn,m 


0.3 


0.4 


0.5 


0.7 


0.8 


0.5 


0.6 


0.8 


0.9 


Case 
(i) 


(a): (>*,«*,) 


(0.001,0.001) 


(0.001,1) 


(0.01,1.999) 


(1,0.001) 


(1,1) 


(1,1.999) 


(1.999,0.001) 


(1.999,1) 


(1.999,1.999) 


0>) : «.Wm) 


(0.1,0.1) 


(0.1,1) 


(0.1,1.9) 


(1,0.1) 


(1,1) 


(1,1.9) 


(1.9,0.1) 


(1.9,1) 


(1.9,1.9) 


Case 
(ii) 


( a ) : z n,m 


0.301 


0.501 


0.601 


0.699 


0.401 


0.599 


0.801 


0.501 


0.601 


(b): < m 


0.4 


0.4 


0.7 


0.6 


0.3 


0.7 


0.9 


0.4 


0.5 


Case 

(in) 


( a ) : C,m 


0.299 


0.401 


0.499 


0.701 


0.801 


0.501 


0.601 


0.801 


0.9 


(b): i*, m 


0.4 


0.5 


0.4 


0.6 


0.9 


0.4 


0.5 


0.9 


0.8 



Table 2: Error due to perturbation 



Perturbation 


Maximum Manhattan Metric 


Error 


Case(i) 


(a): Perturbed (^Jl, !/ m ) 


0.002 


0.0217 


(b): Perturbed (x n ,y m ) 


0.2 


2.1667 


Case(ii) 


(a): Perturbed z n m 


0.001 


0.0227 


(b): Perturbed z* 


0.1 


2.2667 


Case(iii) 


(a): Perturbed t n m 


0.001 


0.0213 


(b): Perturbed ** m 


0.1 


2.1333 




Figure 1: Original Surface 




(a) Simulated Surface due to (b) Simulated Surface due to 

perturbation in independent perturbation in independent 

variables (c.f. Case (i(a))) variables (c./. Case (i(b))) 



Figure 2: Simulated Surfaces due to perturbation in independent variables {c.f. Case (i)) 
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(a) Simulated Surface due to (b) Simulated Surface due to 

perturbation in dependent vari- perturbation in dependent vari- 

able (c./. Case (ii(a))) able (c.f. Case (ii(b))) 



Figure 3: Simulated Surfaces due to perturbation in dependent variable {c.f. Case (ii)) 




(a) Simulated Surface due to (b) Simulated Surface due to 

perturbation in hidden variable perturbation in hidden variable 

(c.f. Case (iii(a))) (c.f. Case (iii(b))) 



Figure 4: Simulated Surfaces due to perturbation in hidden variable (c.f. Case (in)) 
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(b) Simulated Surface due to 
perturbation in all the variables 
(c.f. Case (i(b)),(ii(b))(iii(b))) 

Figure 5: Simulated Surfaces due to perturbation in all the variables 

6 Conclusion 

The present paper explores the stability of CHFIS when there is a perturbation in indepen- 
dent variables, the dependent variable and the hidden variable. The stability during the 
perturbations in all the variables simultaneously is observed to be the combined individual 
effect of perturbations in each variable on the stability of CHFIS. The bound on error in 
the approximation of the data generating function by CHFIS is described individually for 
each case of perturbation in independent, dependent and hidden variables. These bounds 
together are employed to find the total error bound on CHFIS when there is perturbation 
in all the variables simultaneously. The stability results found here are illustrated through 
a sample surface. Our results are likely to find applications in investigations concerning 
texture of surfaces of naturally occurring objects like surfaces of rocks [10], sea surfaces [9], 
clouds [H] etc. 

Acknowledgement: The author Srijanani thanks CSIR for research grant (No:9/92(417)/2005- 
EMR-I) for the present work. 
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